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Abstract 

A compressible turbulence model based on the k — ( (enstrophy) model of Robinson et 
al. is developed. The model incorporates a new compressibility correction model that is quite 
different from previous models. Results are compared to the mixing layer data compiled by 
Settles and Dodson. It is shown that the predictions of the theory are in good agreement 
with experiment. Moreover, the present model is superior to the k — e model with traditional 
compressibility corrections. 


Introduction 

Compressible flows require for their description a velocity field and two thermodynamic 
variables such as the density and temperature. Because of this, the fluctuations of these 
thermodynamic variables are as important as those of the velocity in determining the re- 
sulting turbulent flow. As a result, traditional two-equation and stress models 1 have proven 
to be inadequate in describing such flows. Therefore, it appears that an appropriate com- 
pressible turbulent flow model of the “two-equation” variety, should include six equations 
that describe variances of velocity, density and temperature and their respective dissipation 
rates. 

A simplification of the above approach can be achieved by invoking Morkovin’s hypothesis 2 . 
According to this hypothesis, the pressure and total temperature fluctuations are small for 
non-hypersonic compressible boundary layers with conventional rates of heat transfer, i.e. 

? p«l , |«1 (1) 
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As a result. 
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In the above equation. p'.T^p' and T', u' are the fluctuating density, total temperature, 
pressure, temperature and velocity, M is the Mach number and 7 is the ratio of specific 
heats. The remaining variables represent the mean properties. Based on Eq. (2), equations 
governing variances of turbulent quantities can be taken as the equation for the turbulent 
kinetic energy. However, equations governing the dissipation rates of the resulting variances 
may not be the same. 

The present work is based on an extension of the k — £ model of Robinson et al . 3 to 
compressible turbulent flows. It is shown in Ref. [3] that the k — ( model reproduces all 
available measurements of growth rates and turbulent shear stress distributions for a variety 
of free shear layers, and as such, it should provide a good basis for the current model. The 
compressible enstrophy equation can be written as 
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Here fi; and u- are the mean and fluctuating vorticity, 77 ,- and r- are the laminar and 
turbulent (Reynolds) stress and 6ij k is the permutation tensor. 

The Favre-averaged turbulent kinetic energy equation can be written as 4 
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where Uj and r tJ - are the laminar and turbulent (Reynolds) stress. The dissipation rate is 
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The second term in Eq. (8) can be written as 
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and is normally neglected, consistent with the assumptions of homogeneous turbulence. In 
this work, the term is included in the diffusion term. 

The approach employed in modeling the various terms follows closely that employed in 
Ref. [3] and will be detailed in the final paper. Moreover, model constants developed in 
Ref. [3] were unchanged in this model. 

The compressibility term in the Adequation has generated a great deal of discussion lately 
and led to the developments of models by Sarkar et al. 5 , Zeman 6 , Wilcox 4 and Ristorcelli 7 
among others. None of these models were adopted here because they do not blend smoothly 
with the incompressible limit. Moreover, there is no reason to assume that the time scale 
governing this term is dependent on (. Instead the term 

D \p (<.) 2 ( 10 ) 

is modeled as 

C- (11) 

t p 

where C is a model constant and t p is a time scale. Because the term under consideration 
is zero when the density is constant, t p is modeled as 


— = T\/( Vp- Vp) k 


( 12 ) 


Results and Discussion 

To validate the present model, comparisons were made with the two sets of mixing layer 
experiments that survived the scrutiny of Settles and Dodson 8 . The first set is due to Goebel 
and Dutton 9 while the second set is that of Samimy and Elliott 10 and Elliott and Samimy 11 . 
Rather than restrict our comparisons to growth rates, we opted to compare the predictions 
of the theory with the measured velocity, turbulent kinetic energy and turbulent shear stress. 
Moreover, comparisons were made with a, k — e model using the compressibility correction 
model proposed by Wilcox. It may be recalled that this model incorporates the best feature 
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of both Sarkar et al. and Zeman’s models. Table 1 summarizes the cases given in Ref. [8] 
and described in Refs. [9] - [11]. 

Comparison of present theory with the k— e model and experiment are shown in Figures 1 - 
7. The computational procedure is similar to that employed in Ref. [3]. In presenting, the 
data, we followed the suggestion of Ref. [10] and scaled y with 5 W , the vorticity thickness, 
defined as 

, At/ = U 1 -U 2 (13) 

9 y/max 

where subscripts 1 and 2 represent the fast and slow streams, respectively. Each figure 
consists of three plots for velocity, turbulent kinetic energy and turbulent shear stress, r ry . 
It is seen from the figures that the present theory is in good agreement with cases 2, 3 and 
4 of Ref [9] and all cases of Refs. [10] and [11], For these cases, the k — t model shows poor 
agreement with experiment. The agreement is not good with case 5 of Ref. [9]. The reasons 
for this are not clear at this time. Note that this case is in good agreement with the k — t 
model. Furthermore, applying the compressibility correction of Wilcox to the k — ( model 
yields worse agreement with experiment than the k — e with compressibility corrections. 

Settles and Dodson tabulate and plot (see Table II and Fig. 1 of Ref. [8]) normalized 
mixing layer growth rates versus convective Mach numbers. The data is collected from 
various sources. As may be seen from the data and figure, there is a great deal of scatter 
rendering such tabulations and plots meaningless. The reason for the scatter can be traced 
to a lack of uniformity in calculating and defining growth rates. For most of the available 
data, the growth rate db/dx -=f b/x, where b is the width of the mixing layer. Moreover, there 
is a lack of uniformity in defining the width of the mixing layer. Because of this, judging 
the worth of a compressible turbulence theory based on its prediction of growth rates is 
somewhat misleading. 

The final version of the paper will contain the details of model employed and additional 
comparisons involving wall-bounded flows. 
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Tables 


Table 1: Summary of Modeled Experiments 


Case, Ref. 

u 2 /u x 

P2/P1 

T 0l , T 02 

Mi, M 2 

M c 

2,(9] 

0.57 

1.55 

578, 295 

1.91, 1.36 

0.46 

3, [9] 

0.18 

0.57 

285, 285 

1.96, 0.27 

0.69 

4, [9] 

0.16 

0.60 

360, 290 

2.35, 0.30 

0.86 

5, [9] 

0.16 

1.14 

675, 300 

2.27, 0.38 

0.99 

1, [10] 

0.36 

0.64 

276, 276 

1.80, 0.51 

0.51 

2, [10] 

0.25 

0.58 

276, 276 

1.96, 0.37 

0.64 

3, [11] 

0.25 

0.37 

276, 276 

3.03, 0.45 

0.86 
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Figure 1: Turbulence Model Comparisons with Experiment (Goebel and Dutton Case 2) 
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Figure 2: Turbulence Model Comparisons with Experiment (Goebel and Dutton Case 3) 






Figure 3: Turbulence Model Comparisons with Experiment (Goebel and Dutton Case 4) 
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Figure 4: Turbulence Model Comparisons with Experiment (Goebel and Dutton Case 5) 
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Figure 5: Turbulence Model Comparisons with Experiment (Samimy and Elliott Case 1) 
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Figure 6: Turbulence Model Comparisons with Experiment (Samimy and Elliott Case 2) 
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Figure 7: Turbulence Model Comparisons with Experiment (Elliott and Samimy Case 3) 
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Introduction 

Existing k—u and k—e turbulence closure models are incapable of describing wall bounded 
and free shear flows using the same set of model constants and boundary conditions. It is 
generally agreed that the problem with existing models comes from the dissipation equation, 
e or u. The reason for this can be traced to the fact that these equations were developed 
for high turbulent Reynolds numbers, Re t , but are being employed in situations where Re t 
is low. 

For most turbulent shear flows Re t is typically large enough so that one can assume 
the small scales are nearly independent of the large scales and thus their dissipation rate 
is isotropic. If one makes this isotropic assumption then the terms in the exact dissipation 
equation which depend on the mean flow can be neglected. However, for most shear flows, 
the turbulent Reynolds number varies from very large values in the bulk of the flow to very 
small values near walls and outside of the boundary layer. 

For regions of low Re t , the small scales of turbulence are weakly dependent on the large 
scales, and therefore, on the mean flow. Because of this, terms appearing in the exact 
dissipation equation which depend on the mean flow cannot be neglected. 

The above considerations were the basis of a new two-equation model l . Instead of 
modeling the exact equation for dissipation, attention was focused on the enstrophy or the 
variance of vorticity equation 2 . The model developed in Ref. [1] was used to describe free 
shear layers (wakes, jets, and mixing layers). The same set of model constants was used in all 
calculations. Excellent predictions of growth rates, shear stress, and velocity distributions 
were obtained. 

The object of this work is to use the same set of model constants developed in Ref. [1] 
to study wall bounded flows. The implementation of the model will proceed in a number 
of stages. In the first, the model is implemented in a boundary layer code 3 and the results 
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are illustrated by calculating a flow past a flat plate and its wake. Second, the model is 
implemented in a Navier-Stokes code 4 and the results are illustrated by calculating flows 
past airfoils. The results presented in the abstract are intended to demonstrate that one set 
of model constants can be used to describe all types of shear flows. More results will be 
given in the Anal paper especially those involving adverse pressure gradients. 


Analysis 


The turbulent kinetic energy and the exact enstrophy or variance of the vorticity equa- 
tions can be written as 
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p is the density, p is the viscosity, u\ and u>[ are the fluctuating velocity and vorticity and {/,- 
is the mean velocity. Equations ( 1) and ( 2) were modeled in Ref. [1] and the results can 
be written as 
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with p t being the eddy viscosity. The model constants are tabulated in Table 1. 

Two departures from earlier modeling should be noted. The free shear flows considered 
in Ref. [1] were all characterized by a constant pressure. In the proposed applications the 
pressure gradient is not zero. The pressure fluctuations appear explicitly in the k but not 
the ( equations. Therefore, the turbulent diffusion term in the ^-equation is modeled as 
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and ex p is a model constant. 
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The second adjustment concerns the dissipation term in the ( equation. The turbulent 
time scale ^ can not be less than Kolmogorov’s time scale. To allow for this result, the 
dissipation term now has the form 


For the results presented here 6 = 0.1. 
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Results 


For high turbulent Reynolds numbers, the eddy viscosity is chosen as 


= 0.09 


For bounded flows the eddy viscosity and the turbulent Reynolds number are both zero 
at the wall. Because of this, the above expression is traditionally multiplied by the function 
Jut. There is no unanimity on the form of f ti in the near wall region 5 . More recent imple- 
mentations such as that of Speziale et al 6 employ a function that depends on both Ret and 
y + where 
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where t w is the wall shearing stress. Because separation may take place in the presence of 
adverse pressure gradients a different representation that does not require y + is employed. 
The resulting expression can be expressed as 

/ M = min(/;, 1.0) . fl = ^1 + j exp ( 13 ) 

and 
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The boundary conditions for the two-equation model are 


dC 

dy 


k w = 0 

k„ = ^(TU„f 


U = o 

d 

or — 

(ldk\ 


dy 

\3dy) 


= c 


at 


y = 0 


(15) 


where T is the turbulent intensity (« 1%) and (oo is determined by specifying a free stream 
ratio of ~ 1^ along with k^. For the current work we are using a first order difference for 
|u, which amounts to simply extrapolating £ to the wall. 

All of the results presented here use the first boundary condition for (. Comparisons 
with the other boundary condition will be given in the final paper. 

There are a number of tests a turbulence model must meet to be considered a successful 
model. First, it must predict the correct friction and pressure coefficients and other near 
wall measurements 5 of shearing stress, kinetic energy, and dissipation. Second, it must 
predict the correct growth rate, asymptotic shearing stress distribution, and velocity in the 
far wake. Finally, it must predict the correct ’B’ constant that appears in the expression for 
the velocity in the log-law region, and the manner in which k varies with y in the near wall 
region [See Table 4.5 in Ref. [7]]. 

Figures 1- 3 show plots of k + , ( + , and (— pu'v') + vs. y + in the near wall region where 
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Also shown is the "average” of available data 5 together with the measurements of Schubauer 8 
and Laufer 9 . As seen from these figures, the predictions compare well with the results sum- 
marized in Ref. [5] and experiments. It appears from Figure 1 that the peak k + is under pre- 
dicted, nevertheless, it lies within the range of available experimental data (£ + ss 2.8 — 6.0). 
Figure 4 compares the calculated skin friction coefficient with the correlations of Cole and 
Schoenherr 10 . Again good agreement is indicated. 

The constant ’B’ in the log-law correlation 


u + = —lny + + B 
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is typically 5. Current calculations show the : B ! is not exactly constant, but rather, varies 
slightly from approximately 4.7 to 5.2 over a y + ranging from 50 to 500. 

The next set of figures compares calculated and measured wake data. Figures 5- 7 show 
a comparison of numerical, asymptotic, and experimental results for the defect velocity at 
the centerline, ^ [ U e is the edge velocity, W is the defect velocity], wake half width , b, and 

the maximum where r is the shear stress, as a function of | where 6 is the momentum 

thickness. The asymptotic solution is that of the far wake with a constant eddy viscosity. 
The experimental results are obtained from Pot 11 and Weygandt and Mehta 12 . Again, good 
agreement is indicated. 

Figures 8- 9 compare the defect velocity, W, and the shear stress distributions across the 
wake. Here, also good agreement with experiment 11- 13 is indicated. 

The remaining figures show the Navier-Stokes solution for a NACA 0012 airfoil. Fig- 
ure 10 shows the grid (239x101) which has an initial normal spacing of 2.e _6 C and an outer 
boundary of 15 C. Figures 11- 12 show the pressure and skin friction coefficients compared 
with the experimental data 14 and with the standard k — u> model of Wilcox ' . Once again, 
good agreement with experiment is shown. 


Concluding Remarks 

From the results indicated here, it appears that the model developed in Ref. [1] can 
be used for wall bounded shear flows and it has also been shown that the model can be 
incorporated into either boundary layer or Navier-Stokes codes. The next task is to test the 
model for a variety of other flows including separated flows. 
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Table 1 . Enstrophy Equation Model Constants 
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Figure 1 . Near Wall Behavior of Turbulent Kinetic Energy (Flat Plate) 



Figure 2. Near Wall Behavior of Enstrophy (Flat Plate) 
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Figure 3. Near Wall Behavior of Reynolds Stress (Flat Plate) 
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Figure 4. Skin Friction Coefficient (Flat Plate) 
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Figure 5. Streamwise Variation of Maximum Velocity Defect 
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Figure 6. Streamwise Variation of Wake Half Width 
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Figure 7. Streamwise Variation of Maximum Shear Stress 



Figure 8. Wake Velocity- Defect Profiles 
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Figure 9. Wake Shear Stress Profiles 



Figure 10. Partial View of C-mesh for NACA 0012 Airfoil 
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Figure 11. Pressure Distribution for NACA 0012 Airfoil (M 0 0 = .502,i?e oo = 2.91e6, a 
2.06) 



Figure 12. Skin Friction Comparison witb k to for NACA 0012 Airfoil (A^co — .502, Hcqq 
2.91e6, a = 2.06) 
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